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A STUDY OF THE DISTRIBUTION OF MOLECULES 
UNDER FREE MOLECULAR FLOW CONDITIONS 
AFTER COLLISIONS WITH SIMPLE GEOMETRIES 

By 

David W. Tarbell and James 0. Ballance 
SUMMARY 

The study makes a detailed examination of molecular distribution 
under free molecular flow conditions within surfaces of simple geometries. 
By considering infinite length in one dimension, the problem could be 
reduced to a two-dimensional analysis. Exact, closed form solutions were 
obtained only for circles; however, numerical solutions were obtained for 
ellipses and parabolas. Monte Carlo computer technique results were 
obtained and compared to the theoretical solutions. Derivations of (a) 
the distribution of molecules over the surface as a function of collision 
number, (b) the exit distribution across the opening of each surface con- 
figuration, and (c) the distribution across the center line of each con- 
figuration are presented. 


SECTION I. INTRODUCTION 

Studies of free molecular flow problems have usually been concerned 
with the gross effects resulting from the flow such as the conductance of 
a tube or the forces on a body moving at a velocity large compared to the 
thermal velocity of the gas molecules. The analytical solutions are 
usually obtained in exact forms for limited geometries. .The use of high 
speed digital computer methods [l, 2, 3, 4] has allowed the determination 
of many interesting parameters such as distribution of molecules, average 
number of collisions within surfaces, etc., parameters which would have 
been difficult if not impossible to determine from exact analytical 
studies. It was felt by the authors that a detailed examination of 
molecules under free molecular flow conditions would allow much greater 
insight into actual physical processes and would assist in practical 
problems such as cryogenic pump array design, molecule -surface inter- 
action studies, high altitude atmospheric measuring programs, etc. 



Typical of most studies of free molecular flow, two assumptions have 
been made: (1) The molecules interact only with the surfaces and (2) 

molecules colliding with the wall are re-emitted diffusely, (The terms 
"molecules" and "particles" are used interchangeably throughout.) In 
all cases, the surfaces examined are considered infinitely long in the 
"z M direction such that the analysis can be reduced to two dimensions.- 
Theoretical studies presented in this paper consider only elementary 
surfaces which, in two dimensions, are circles, ellipses, and parabolas; 
however, completely digital computer studies have considered hyperbolas 
and triangles [3]. Only for the case of the circle could solutions be 
made exactly in closed form. The other cases resulted in integrals, 
some of which were solved by computer methods. Monte Carlo methods used 
for this study are described in Reference 3. 

This paper presents the derivation of (a) the distribution of 
molecules over the surface as a function of collision number, (b) the 
exit distribution across the opening of each surface configuration, and 
(c) the distribution across the center line of each configuration. The 
results of numerically integrating these expressions are given and com- 
pared with the results obtained by the Monte Carlo method. 


SECTION II. DISTRIBUTION OF MOLECULES OVER A SEMICIRCULAR SURFACE 

AS A FUNCTION OF COLLISION NUMBER 



FIGURE 1 


3 


c 


Consider .a flux of molecules entering the semicircular geometry 
shown in Figure 1. After being reflected at the surface, some of the 
molecules will exit the array, while others will collide again with the 
surface a number of times before exiting. 

The problem to be considered here is to find expressions for the 

distribution of these molecules over the surface, as a function of 

• * 

collision number. When the initial distribution is known, the problem 
is simply a matter of trigonometry and calculus, as follows: 

Let the initial distribution be given by 

= number of molecules per radian that collide with the 

surface at point P on the first collision. (See Figure 1.) 

Then the total number of incident molecules, N Q , is given by 



o 


According to Lambert’s Cosine Law of Diffuse Reflection, the number 
of molecules per radian that are reflected from Pq at an angle <p with 
the normal is proportional to the cosine of this angle, i.e.. 


%,cp ■ k cos 9. 


Then clearly, 


% = 



cos cp dcp. 


where Nq; is the surface distribution function for any arbitrary collision. 
This yields immediately 


t 


N, 


a,<p 


1 

2 


N a cos 


cp. 


so that 



4 


Now the number of molecules that collide again with a section of the 
surface between (3 and rt/2 (shaded band in Figure 1) will be the sum of 
three contributions: (1) for OL > it/2, (2) for it/2 > O, > p , and (3) 
for P > a. 

Contribution (1): Ot > rt/2 

The number of molecules per radian reflected from Pq, that hit 
between P and it/ 2 is 


2% cos cp dcp 


1 __ r ■ 1 e + 

2 N a [sxn cp] ^ 



5 


so that 






Then the total number 
from all Pq, with OC > jt/2. 


of molecules that hit between p and tt/ 2, 
is 




it 



Jt 


2 



da . 


(l) 


Contribution (2): 2 — ^ ^ 

In this case there will be two contributions: (2-1), angles greater 

than a, and (2-2), angles less than a (Fig. 2). For the first contri- 
bution (2-1): 



Figure 2 


sin (-6) 


N a,p = j 2 Na cos 9 dcp = 2 N « 

« 

" 2 

N a"p = \ % O " sin 5 >> 


but, from Figure 2: 



so that 


2 “ 1 1 

N a,p “ 2 Nq 


1 - sin (f + §) 


For the second contribution: 


it/ 2 


2-2 P i l / 

N ry R ' o N a cos qp dcp = - N a ( sin -r - 


7 


or 



Therefore, 


p 2-i 2-2 i 

N 2 = M +N = — N 

w a,p .a,p a,p 2 a 



Then the total number of molecules that hit between p and it/ 2 9 from all 
Pq, such that it/2 > Ot > p is 



Jt/2 

j | n q [2 - sm (f + f) 
P 


COS 






da 


Contribution (3): 0 — P 


From Figure 3: 
e+& 

” J I N a cos ^ d<p ’ 

e 

N a,e ■ I “a < sln <e+5) • sln e) - 


( 2 ) 



8 



Now 


2(e4€>) - * - (p-a) 


€ + 8 = f - (^r) 


sin (€+6) = cos 


p-a 


Also 



so that 


9 



and the total number of molecules that hit between (3 and tt/ 2, from 
all Pa, such that f3 > a, is 



o 


Combining equations (1), (2) and (3), the total contribution to 
the number of molecules hitting between (3 and n/2 from all Pa> is: 


Ti 



(3) (A) 



o 


( 5 ) ( 6 ) 

Now terms (2) and (4), and terms (3) and (6) can be combined directly. 
Also, since 




the expression reduces to 



P o o 


Now, from the symmetry of the problem (Figure 1), it is clear 

that 


N *-a = N a 


Also 



and 


d(rt - a) = - da, 

so that by substituting (n - a) for a in the first integral: 

/ \ K a sln (f + f) da ■ / I \-a sln (f + <■< * - a) 

:/2 n/2 

Jt/2 

= f sin (t + da. 
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so that it cancels with the fourth integral. With this simplification, 
the expression becomes 


it rt/2 

B P - -/ \ “a “ s (^r 6 ) d “ + /»« da + 

P P 



(4) 


This is the total number of molecules that hit between p and ir/2*from 
all points on the surface. Since this will vary with the collision 
number, n, let NX n ) = total number of molecules that hit the surface 
between (3 and rt/2 on the n collision. 

Now the number of molecules per radian that hit a point , P^,, on the 
n fc h collision is found by calculating the limit: 1 


lim 

A|3 -4 


N 


(n) 


- N (n) 


Ap 



But this is clearly just the derivative 


-d 

dp 


; therefore, 

P 


N 


( n ) = zl N ( n ) 

a dp p 


p = a 


so that 


*/2 


n' n ' 1> da - 


~-«r] 


it/ 2 


From Figure 1 it is clear that N^ 11 ^ = 0 for p = rt/2 and any n 


Therefore , 


it/2 



P 


da = N 


(n-l) 

P 


(5) 
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Also , 


cos 


a - p 


g 6 , . a . 6 

= cos — cos 2 + sm — sm ^ 


( 6 ) 


Now, let 


F (n) (a, P) = \ 


_ 7 (n-i) ( a 6 , . a . b\ J/v 

' [cos — cos ^ + sin — sm ^ ) da. 


(7) 


Then, using equations (5), (6), and (7) in equation (4), it follows 
that 

N ( n) = . F (n) (it, p) + F (n) (P, P) +N^ n " 1) + F (n) (p, p) - F (n) (0 , P), 


or 


N< n) - 2F <n) 0, P) - F (n) (0, P) - F <n >(«, P) + N*”- 1 ) (8) 


and 

■ (9> 

p « a 

which is the final expression for the distribution. Given any initial 
distribution, N^) , the distribution at the next collision can be found 
by first calculating f(2) (OL, p) from equation (7), using this expression 
(with the indicated substitutions) in equation (8), and then using 
equation (9). 

For specific cases, the elementary integrals given in Appendix A 
will be useful. 


13 


B. Case I: Incident Distribution Uniform over the Surface 

1. Distribution at Second Collision 

For an incident distribution of molecules uniform over the 


surface , 


/IX N 

(1 ) o 

Nj, = — * constant, 
CC xt 


( 10 ) 


where N 0 is the total number of incident molecules, 
equation (5), 


Then, from 


it/ 2 


N 


(i) 


-I 


Till 


N a * da 


N 


da. 


„(1) _ N 0 fit o , 


(ID 


From equations (7) and (10), 

a 


F (2) (a, p) = \ 


p No/ a 3 , . a . b\ _ 

/ -£• f cos cos “■ + sin j sin ^ J da, 


F (2) 

(a, 

P) 

N 

o 

~ 2jc 

Then 

* 



f (2) 

(P, 

P) 

N 

o 

2jc 

f< 2 > 

(0, 

P) 

N 

o 

2tc 

f< 2 > 

(*> 

P) 

N 

o 

2jt 


a B „ a . B' 

— cos ^ “ 2 cos 2 Sin 2 


,2 sin | cos f ~ 2 cos \ sin \) ~ 


N o . £ 

T sin v 


%- 

= L,e 


2 ’ 



and, using equations (8) and (11): 


b < 2 > . 2F (2) (P, P) - f < 2> (0, p) - F (2) (n, p) + n* 


N o . p N o e/.f'i, A 

= — sm 2 ■ ~ cos o + T l “ " P )• 


2*^2 


N p 2) = T (f • p + sin I • cos f) 


From equation (9): 




-1 + J cos ^ + 


1 sin fj > 


N o /. a . a 

- u{ 2 - COS 2 + Sln 2 


2. Third Collision 


From equations (7) and (13): 


f ( 3 V, p) 


„ a . a\ / a 6, .a. 

2 - cos — ” sin j ( cos — cos + sin j sin 


F (3) (a, p) 


o , . a 6 . a . 8 /a , 1 . .A 

•j— 4 sm — cos 2 “ ^ cos 2 sin 2 \2 + ~2 sin y cos 


- sin 2 ^ sin t _ s in 2 ? cos ^ - (% - 4 sin sin tl. 



15 


Then, 


F ( 3 ) ((3, p) 


= £ [" (I + \ Sin P ) C0S I 

- (f ~ \ sin.p^ sin |j 


, 2 P • 6 

sin f sm f - 


sm ^ 


N r 
o 

4jt 


- ^ + j sin p) cos ^ - sin ^ + cos 2 ^ sin ^ 


. 2 B 6 / 8 1 . Q \ . B _ 

in £ cos f " ^2 ~ 2 Sin ^ ) sin 2 


sin 


4* L \2 2 


sin pj cos ^ - sin ^ ~ cos ^ sin p 


1 . 6 . Q 
2 sin 2 sin P 


(f ' \ sin P ) sin I * 


N r 


f( 3 ) (P, P) = 4 ^ 


2 (cos & + sin - sin ^ 


F (3) (0, P) 


= £ (-4 sin |). 


F (3) (it, p) 


N ° ( . 6 it 6.6 6 « . 6’ 

4 ^ 14 cos f “ J cos 2 ’ sin 2 “ COS 2 “ 2 Sin 2 


Then, using equations ( 8 ) and (12): 


n< 3) - 2F (3) (6, 6) - f ( 3) (0, p) - f ( 3) (it, p) +N^ 2) 


4jt P V C ° S 2 + Sin 2 ) 


- 2 sin ^ + 4 sin 2 " ^ cos 2 


4* ^ f cos *| + sin + cos & + sin ^ + 4 sin ^ - 4 cos ^ + 2xc - 4(3 


N p 3) = 2^ _ 2,t - 4 e + (f ■ p + 7 ) sin f + (f - p - 7 ) cos I 


From equation (9): 


,(3) = . JL n O) 

a dp p 


* 4 + Kf ■ p + 7 ) COS f ■ Kf ■ 


p - 7 J sin ^ - sin ^ - cos | 


4 3> = 5? 8 ' (f • “ + 5 ) cos f + (f - a - 5) sin f . 


3. Fourth Collision 


From equations. (7) and (15): 


r <4> («, »-& 


8 - ^ - a + 5^ cos j - a - 5j sin 


( a £ 

cos — cos . 


+ sin sin ^ da 
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IS [l 6 (sin § cos 


£ - cos J sin f) " (f + (f + 7 sin °D cos f 


^ + 5 ^ sin 2 j sin ^ ^ + a sin a + cos cos ^ 


+ f sin a - a cos Cy sin ^ + f j - 5 y sin 2 j eos 


f + (f ‘ 5 ) (f ’ i sin “)• 


sin^ - j (^ sin a " a cos °^) cos 2 ~ 2 (^T~ ~ 


a sin a - cos a sin 


,ln £ 


Then 


r r \ N 

F< ’ <». » - is 


‘ (f + 5 ) (f + 5 8ltl B ) cos f ' (f + 5 ) sin2 f sln 


2 £ sin ^ 


i ref. 


+ | l ' + P sin p + cos p 


S ) COS 2 + \ ( 


sin p - p cos p sin 


£ 


+ lf - 5) sin 2 | cos f + (§ 


5 ) (f ” 2" sin sin 2 “ J ( Sin P " P cos p) cos 2 


- 2 {% P sin P ~ cos Pj sin 2 


N 

o 

16 * 


cos 


I [- f (2 + 5 ) + \\T + p sin p 


+ COS P I - 


- [ sin p - p cos P 


^ + sin jT- (j^ + 5^ + ^ ^sin p - p cos p^ 



5 ) | - y My- - P sin p - cos p 


1 


•|cos f f (jr + 5^) + f ^1“ + 2p sin f cos f + 1 - 2 sin 2 |^) 


2 sin j cos 


| - P + 2p sin 2 0j + sin f 

in ^ rns ^ « 9ft nne^ ^ 


-if + 5 


2 sin | cos f - 2|3cos 2 * + p f ' 5 ) f 



^ — 2P sin f cos f - 2 cos 2 f + 1 


)]} 


COS 


e 


-!(! + 5 ) + Kf + 1 ) + f 


+ sin f 


f + 5 ) + f 


Of-Kf *■)]} 


P, P) = 


N r 

= -2-1 

r. i/j + 4 

-§) 

+ il 

cos ^ + 

16rt 

2 \2 

V 

2 



£ (l - 4 - ^ - * -ill 

2 V7 2/ 2 2 


i\ rt 111 . 6\ 
'.) " 2 2J Sin 2j 


°> P > = I§^ (- 16 sin f + f cos f + f sin f) 


r 16 cos f - + 5^ f$\os | -(x ■+ 5^) sin f f 


- 1 ) cos ^ 


*, P) 
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i <*> sl ” f + (f - 5 ) cos I + (f ■ 5 ~) (f) sln f - 1 <*> f 



s in & 

V 2 ) 



f(4)( *- p > = if; [( 10 § - ¥) “• I + (- 4 • ^f) sin fl • 


Then, using equations (8). and (14): 
(4) 


Np ' - 2F (4) (P, p) - F (4) (0, p) 


f ( 4) (*, p) + n< 3 > 


N 

o 

16* 


+ f- + l 


¥) 


\ • 10 I + T^ COS 2 


§ 


+ - 4p - - jt - 11 + 15± + 5^ + sin | + 


4^| - p + 7^ sin ^ + 4^| - p - 7^ cos ^ 


8« - 16p 


N, 


(4) 


' T§7 [ 8 ” ‘ 16p + (’ 38 + T ‘ 8p ' ? + ?) cos I 


+ (38+^- 8 <s+f -£).!„§ 


From equation (9): 
,<*> = - 


N. 


a 


± N wl 

dp N J 


p=a 


( 16 ) 
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- 16 +(- 8 - § + p) cos f + (- 38 + y 


2 2 


+ 8 + sin ^ ^ ^38 + y - 8|3 + 


& . §l\ 

2 2 ) 


4 4) - 3^ 32 - r. 22 + 5 f - 6« + 


ita af\ 

2 " 2 / 


a 

cos 2 - 


L 7 it , , jta cA . a 

-[22 - — + 6a + -^ 2 ~J sin — . 


4. Fifth Collision 


From equations (7) and (17): 


F<5, <“, » - m 


32 - f 22 + y - 6a + y - t) cos f 


- ^22 - ^ + 6a 4- -y - Yj sin J c°s “ cos ^ + sin j sin -| j da. 


F<5><a - P >’6fel COS 


| [64 .inf -(22 + $(§ + i sin J) 


+ ^6 - y| fj) fy + a sin a + cos qN ) + 7 (^y + (y - 1 ) sin a + a cos a) 


2 V 6 ' v 2 


+ 0r - 22^ sin 2 J ■ (^ 6 + f) (j) ^ sin a " a COS 
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1 l ~ 

+ — (OL sin OC + cos CL - ^ cr cos a) 


+ sin | 


- 64 cos ^ - [ 22 + ~ ) sin 


+ ( 6 - (sin a - a cos a) + i sin a + cos a - ~ cos 


a) 


(l— 

V 2 




sin a - cos a 


•) 




1) sin Q! - CC cos a 


Then, 


,(5) 


N 




<P - P) " 6 S 1 cos 2 


64 sin 


£ . 

2 


22 


*¥ 


£ 

2 


sin ^ 


cos 


I 


D 


+ 6 - i 


V w v 


+ J + (^- - 1) 2 sin | cos -| + p - 2p sin 2 ^ - 22^ sin 2 £ 

- ^6 + (jJ ^2 sin ^ cos £ - p + 2p sin 2 ^ + 7 (^ 2 P sin ’f cos 7 + ^ 


- 2 sin' 


2 £ _ §1 

2 2 


+ p 2 sin" 


£ 


+ sin | 




( 6 ■ -fX*X 2 sin 


e 

2 


cos 


£ 

2 


2p cos 2 | 



N>IQ 






23 



N>|S* 



- - 152 + 18* - 32(3 - 2itp + 2p^ + sin ^ 51 - — - I5p 


3 *P - ^ f + 60j - f + 25 ±- + lo| * - £ + §5 + 152 + 14 * 


- 32(3 + 2rtP - 2(3 J + 32rt - 64p 


N p 5) = ^ [ 32 * ’ 64P + (- 187 + 28 J« - V + " ^P - 4*(3 + ^ 


- ^ + %r) cos | + f 187 + 18^ tc - ~ + £ - 47p + 5itp - 2|1 - ^ 


4 6 


2 _ q 2 


2 4 24. 


2 4 


+ SlA sin | 


From equation (9): 


- — H 

ap "p j 


6^ ‘ 64 + (- 2 


187 + 28| jt - + | 


- 47P 


4jtp + 2|f .«|f + ^) sin | + (. 47 


_ j£ + £l 


4it + 9p - * + ^ ) c 


4 fl87 + 18^t - + €: - 47p + 5rtp - ^ ^ cos f 




N 

o 

128* 


128 + 


1 * 2 

93 - 10 — * + 



29a - 4*a + 


7a 2 


Tta 2 a£\ 
4 ’ 6 J 



93 


+ i 4 



29a 


3*a 



jta 2 

4 




(19) 


C. Case II: Incident Distribution Uniform across the Y-Axis 

1. Distribution at Second Collision 

For an incident distribution of molecules uniform across 

the y-axis: 


„ dN , . . 

N a = d5 4 constant 


but 


Now 


— = constant = k. 
dy 


_ dN _ dN # dy 
N a “ da ” dy * da * 


and from Figure 1 it is clear that y = - cos a so that 



Therefore, Nq; = k sin OL , but 


N 

o 


/ 


k sin a da = 2k 


so that 


N a = 2 N o sin a • 


From equation (5). 


Jt/2 


it/2 


,(D 


P 


= da = J |n q sin a da 


p 


= | N o (- cos f + cos (3) 


N p 1} = 2 N o C ° S P 


From equations (7) and (20): 


f ( 2) (a, p) 


i N o 


sin a 


a _ , _ 

cos — cos ^ + sin 


£ 


/r. 


6 

cos £ 


. a 2 

sin — cos 


T + sin | 


sm‘ 


MQ 
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Then 

F (2) (p, P) = y (sin 4 § - cos 4 0 

= y (sin 2 ^ - cos 2 (sin 2 ^ + cos 2 

r?\ N 

f u \p, P) = y (- cos p) 

F (2) (0, *P) = y (- cos f) 

F (2) (rt, p) = y (sin |). 

Then, using equations (8) and (21): 


N< 2) = 2F (2) (P, P) - F (2) (0, P) - F (2) (rt, p) + Np 1} , 
= y (- 2 cos p + cos - sin ^ ^ cos p^, 

„(2) N o ( B .PI A 

N p = — (cos ^ - sin - — cos PJ. 


From equation (9): 


n< 2 > = - JL n (2) 1 
a dp p J 

p=a 


N 

o 

3 


1 . B 1 p , 1 . 0 

2 sin £ " "2 cos 2 + 2 sin ^ 


p=a 


..(2) N o ( , a . a . A 
N a = ~ I sin — + cos — - sin al. 


( 22 ) 


( 23 ) 
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2. Third Collision 

From equations (7) and (23): 

F^(a, P) = yy f fco s — + sin ^ - sin a 


a e . . a . 

cos — cos ^ + sin — sin 


N 

o 

12 


cos 


B /a . i . ^ , .2a, 4 3 a' 

^ - sm a + sin - + - cos J - 


2 ^2 2 


2 3 


. B f . 2 a , a 1 . 4 . 3 a 

+ sin I sm 2 + 2 ~ 2 sin a " 3 sin 2 


Then 


F (3) (f3, P) 


!o r 

12 


COS 


I (^| + j sin P + sin 2 | + j cos 3 


+ sin - fsin 2 \ + \ ~ \ sin P " J sin3 f 


N 

o 

12 


cos 

'6 . 6 
^ + sin *7 cos 

§ + 

* 2 P , 4 
sin *7 + — cos 

§ . 

4 6 

— cos £ 

sin 2 f) 

2 ' 

^2 2 

2 

2 3 

2 

3 2 

2 / 


, . 6/. 2 BJ .6 P 4.6_._4. 6 2 B 

+ sin | 11 - cos 2 + 2 ” Sin 2 COS 2 “ 3 sin 2 + 3 Sin 2 COS 2 


,,, N 

f ( 3 ) (P, P) = 


J cos P + 2 cos 2 + ^ + 2^ sin 2 


£ 

7 


f ( 3) (0, p) 


!° r 

12 


COS 


£ • A 

2 { 3 J 


f (3) (*, P) 


V 

12 


cos f (f + l) + sin I 


l + f 
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From equation (22); 


Np 2 ^ = yj ^4 cos ^ - 4 sin ^ - 2 cos p^. 


Then, using equation (8): 


N< 3) = 2F (3) o, p) - f ( 3) (0 , P) - F (3) (it, p) + n< 2) 


12 [3 


£ 


cos p + p cos ^ + (2 + p) sin ^ J cos ^ 


(f + 1) cos -| 


-(H) 


* A . p , , 6 . . ft 

— - t sm * + 4 cos - 4 sm t - 


2 ~ 2 cos p 


N r 


n ( 3) = 
p 12 


J cos P + [ P 


-&+£) 

2 + 3y 


cos -| + ( p 


-s-£\ 
2 3/ 


sm 


£ 


(24) 


From equation (9): 


N< 3 > = - i. n (3) ] 

a dp lN p J 

p=a 


N 

o 

12 


- sin p + p 


•f + f)(- sln f)(?) 


+ cos ^ + fp - If - f cos sin ^ 



27 \2 


p=a 


“ " if [* f sin p + (' I + t + i) sin f + (f ‘ f + ?) cos f 


"p=a 


N 


N 


(3) = % 


a 


o_ 

12 


( 1 * . a\ a ( 1 it a\ a , 2 . 

6 " 4 + ij Sin 2 + V" 6 + 4 ' 2j cos 2 + 3 Sin a 


(25) 



3. Fourth Collision 


vo|oo 
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f ( 4) (P, P) 







sin p 


+ 


1 

4 



P + 



9 cos 


3 £ 

2 


+ sin 


p 17 i 
2 |\“ 6 




+ 




sin p + 



cos p 


+ 


12 






+ < " \ " f) cos | - 


" 1 " f) COS 2 sin2 f + 


± + z)b 
12 8/ P 8 


COS 




8 26 . 2 B 

- cos ^ sm ^ 


1 , Jt\ . ^ 

' 6 + 4j Sin f 


1 


- + & 
6 ' 


f) sin f cos2 f 


+ - 


sin 2 ^ cos H + 2 


! + £> 
4 k 


. B 2 P 
sm ^ co s ^ j 


-i + fjslnf 


i- - p + EL 

12 87 p 8 



F (4) ((3, P) 


= _£ r cos £ f . £ + *£ . I . £T) + sin £ f i- + * - £ . *B + £- 

COS ~ l /. O I ' OJ - L1 O l 1 O ~ 7. o o ’ o 


2 V 3 8 4 8 


2 \12 4 3 8 8 


- 9 cos p . 


f ( 4 ) ( 0 , p) = 2 I cos f (- \ ' f) + sin f (- £) * 


F ( 4 ) (rt, p) ■ 2^ ] cos f 


i - i) + (- 12 + f) * 


. B 
sin | 


- \ + f) + (j ~ i) + 


1 rt \ if , 8 

12 " 8 / * + 8 + 9 


_(4) , QN o B ( 1 it). . BA. it\ x 8 . 6 

F (rt ’ P) = 24 L COS 2 \J2 - Uj + Sin 2 Vj 2 “ u) + 9 Sin 2 * 


£ rli - ^ + i sin £ 


From equation (24): 




jt + cos *2 + ^ cos P 


Then, from equation ( 8 ): 


n1 4 > = 2F ( 4 ) (P, P) - F ( 4 ) (0, P) - F ( 4 ) (rt, B) + n£ 3) 
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. . P (I , Jt 26 rtB 6~ 1 1 * 

+ sln 2U + 234 + 4 + 4- l2 + 12 


~ - n + 2p 


“ 9 cos p 


^4[cosf(3| + ^ + f-£-ll f ) 

+ sln I (• 4 + ¥ - ¥ + • lf) - 1 c ° s e] 


From equation (9): 


,(*) = 


_d N (4 > 

dp W P 


N r 
_o 

" 24 


gf4 , it 26 

COS 2 3 + 4 4 


- cos ^ T- 3- + ^ ^ ^ 

2 os 2 \ J 9 3 4 4 12 J 


+ ( - I sin fiV 3^ + ^ - ii£ 

V 2 2A 9 + 3 4 4 12 


+ sin | (j- “ f + f) + | sin p 


jt a jta o?\ a . /'ll 5n . a , «a oF\ . a 

24 " 6 + T - “J COS 2 + ll8 • 24 + 6 + T - “J sin 2 


4 

- sin Q 



Fifth Collision 


From equations (7) and (27): 



U rt_ 
18 24 



ita af\ 

8 8 y 


cos 


a 

2 


+ 


11 5* a 

18 24 " 6 


ita af\ 

8 8 ) 



sin a 



0 ! 


2 cos 


o: pa 

+ sin j sin | ) da 


P 

cos J 


11 , Jt \ fa , 1 . 

• 18 + 24^2 + 2 8,1 



*\ /a 2 . a sin a 

8 ) y~ + — — 


+ (rj- - 1) sin a + a cos a\ + y|- + jjf) sin 2 j 


c ^ I sin a - a cos a) + ^fa sin a + cos a - J a 2 cos cc 



• 3 a 

- cos° ■=■ 

* z 


+ sin ^ 


11 . Jt \ . ‘ 2 a . A ic\ A\ 

“ 18 + 2V Sin • 2 + ^6 ■ 8j {?) * 


a cos ^ sin cc + cos a - | a 2 cos 


5jA /a A A jt\ /a 2 a sin a cos a) 

« k-7 ««<»• r + 7 h o o — 


£ 

2 


cos os 
2 , 


35 



36 


2 48 18 24 


1-2 sin 2 + -4 + 4? J sin 2 ^ - «|y cos ^ ( 1 * sin2 ^ 


- & 
*f sm 


■ n t |i\ B + r 1 + £' 


36 48 


24 32 


36 12 8 3 8 


sin •£ cos & + + £& + £■ ) [ 2 cos 2 - l) + 

sm 2 cos 2 V24 24 16 16 16 J V 2 7 48 


+ 1 “ I? + 24 ) 6 - COs2 2j + 17 Sin 2 l 1 ~ COs2 f 



f ( 5) 0, p) = - 


N 


{ cos I 24 “ 16 + (l2 ” 72 ) P " (48 + 32 ) p2 + 48 p3 ] 


+ sin £ 


59_ it /_«_ 7 V TjL it Vg , i_o 3 l 16 / 

72 “ 48 \24 " 72 / \ 48 " 32 / 4 ^ J 27 \ 


16 7 2 6 .26 

" 27 V COS 2 ” Sin 2 


F (5) (0, p) 


N r 
o 

48 | COS 


1 


5_ 

24 


It 

16 


16 

27 


+ sin -| 


I- + iL. 
24 16 


} 


F (5) (n, p) 


■ - it { cos f [ (- 3? + Is) * + (k - ^ 


5. 2L. + jd 11 

24 16 48 18 


+ sin *1 


[7 it + 5*V 

36 48y * 



+ 


(- 




f (5) U P) 




59 35 tc , rtfN 

2 y 72 " 144 16 _ 96j 


+ sin 


£ 

2 


( I7_^j£ jsf\\ 

Vf 216 " 144 ' 16 “ 96 yj 


From equation (26): 


N 


n (4 > = -a 

W f3 48 


COS 


£( 1 ° + 8 

2\9 3 


P + 


2 


£ 


liiA 
6 ; 


+ sin 


£ 

2 


Z£ + 8 jjg 

9' f 3 P 2 + 2 



8 

■9 cos 


P 


Then, from equation (8): 


N^ 5) = 2F (5) (p, p) - F (5) (0, p) - F (5) (*, p) + N< 4) 



5_ 

12 



j_ , it , 16 , 59 , 35 rt 

24 + 16 + 27 + 72 + 144 " 16 + 96 






5_ rt_ 

24 " 16 


, 49 17 n rtf «f 

216 144 16 96 


8 Q rtg . jgf 5 j 

3 P " 2 + 2 6 


. 32 _ 8 

+ 27 COS P - g- CO 


N 

N (5) = — 
W p 48 




A 3 , io3p rt p £f 

T ** A ' A /• » A " A / ■ -I /- “ A I 


16 96 


+ sin 


• 1 f c\l 61 jt , rt 2 

Ln 2 l" 6 TM “72 + 16 


103g 7«g Ugf , rtp 2 P 3> 

' A /■ " 1 A * A / ' 1 /" “ A / 


+ — COS p 


From equation (9): 


,(5) . 


dp p 


_ cos £(l03 «_22g 2se ._3£ 

48 2 36 3 24 16 24 


, 1 6 
+ Y COS 


-17 61 it , rt 2 103§ 7rtg.llgf.rtgf £_ 

T08 '72 16 ' 96 36 12 24 16 "24 


.6 

+ ( *2 ) [ - sin | 


145 it .rtf rtf . 103g «§ llgf .gf . jtgf N 

72 16 " 96 36 3 24 24 16 y 


+ sin | (rn - is + lie + |ae - #3 - |* . ta e 
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n ( 5 ) N n 8 . i jt 2 . it 3 37a _jta 5a 2 ita 2 . o?\ a 

N a = 48l_27^ 1 - na ‘\216 + 144 " 32 + 192 “ 72 + 6 “ 48 " 32 + 48y) COS 2 


, (i 47 61 jr jtf_ 37a jta 5a 2 a 3 itaf\ . a 

l 216 “ 144 32 " 192 72 24 " 48 "48 32 ) n 2 


D. Summary of Expressions for Distribution over a Semicircular 
Surface at the Collision 

1 . General Expr es s ions 


F (n) (a, p) =\ 


x (n-1) ( a 8 . OC . 8\ , 

N v ' cos *r* cos •£ + sm sm aa. 

a , V 2 2 2 2/ 


N ( n > = 2F (n) (p, p) - F (n) (0, p) - F (n) ( rt , p) + N< n_1) 

P P 


,(n) _ 


±- N < n > 
d P w p 


2. Incident Distribution Uniform over the Surface 


f <*-*> 


N e = “ f ? " P - cos ■§ + sin 
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(5) 

a 


N 

o 

128* 


128 + 


93 


21 * *2 

2 T 


24 


+ 29a 


4*a + 



no 2 

4 



• cos § + (- 93 + ^ - ^ + - 29a - 3*a + ^ 

rta 2 a 3 ^ . a 

sin 2/ 


(19) 


3. Incident Distribution Uniform across the Y-Axis 


\ N sin a. 
a 2 o 


( 20 ) 


N 


(1) _ I 


= — N q cos p. 


( 21 ) 


N 


(2) N o / 6 -31 

p = — M COS - sm - — cos e 


( 22 ) 


(2) N o / . a , a 

= — ( sm J + cos -J - sm a 


(23) 


N 


(3) 


N 

o 

12 


(e - 2 ~ j) sin 2 + (& " ! + 4) COS "f + "f cos P (24) 


N 


( 3 > = _£ 


a 


N 

£ 

12 


I . Jtt + a 


, ^ , . a , ( 1 . * a\ a . 2 . 

6 " 4 + 2j sin 2 + r"6 + 4 - l) COS 2 + 3 Sln a 


(25) 
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SECTION III. DISTRIBUTION OF MOLECULES OVER ELLIPTICAL 
AND PARABOLIC SURFACES AT FIRST COLLISION - 
INCIDENT FLUX UNIFORM ACROSS Y-AXIS 

For any given geometry, the number of molecules per radian that 
collide with the surface at an angle a, measured from the negative 
y-axis, is: 


_ dN 
a " da 


dN _ dy 
dy ’ da 


Assuming a uniform incident flux across the y-axis: dN/dy * constant « k, 

where the constant k is easily determined from the condition that a total 
of N q molecules enters the geometry across the opening from y = -1 to 
y = +1. Then 

+1 

N q = J' kdy = 2k, 

-1 

or 


k 2 N o» 


so that: 


N 


a 


1 

2 


N 


d Z 

o da 


(30) 


The derivative dy/da will depend on the geometry of interest. 
A. Ellipse 


Consider now a semi-elliptical surface with its opening, two 
units wide, along the y-axis (Fig. 4). The general equation for such 
an ellipse is: 

X 2 1 

jj + y - i. 


( 31 ) 


FIGURE 4 


(Note that when A = 1, this is the equation for the circle which was 
treated in the previous section). From Figure 4, for a > it/2, 
tan (a - Til 2) = y/-x, and, since tan (a - jt/2) = - cot a , this may be 
written 


y = x cot a. 


For a < Til 2: tan a' = — = — , so that, in both cases; 


-y y 


y = x cot OC 


(3 


and 


x = y tan a. 


(3 


Using equation (33) in equation (31): 


y 2 tan 2 a . „ 2 _ .. 
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or 


tanf _ g 

A 2 


1 


and 


_ ± A ^ _ ± A cot a 

^ (tan 2 a + A 2 ) 1/ 2 (1 + A 2 cot 2 a ) ± / 2 

But, when cot a > 0, y < 0 and when cot a < 0, y > 0, so: 
y = - A cot a(l + A 2 cot 2 a)" 1 / 2 . 


(34) 


Then 

- (- A cot a ) (- '(1 + A 2 cot 2 a)“ 3 / 2 (2A 2 cot a) (- esc 2 a) 

+ A esc 2 a (1 + A 2 cot 2 a) -1 / 2 

: - 1 ). 

y 2 - 1 = - x 2 /A 2 , so that: 


da 

Also 


y esc 2 a f -x 2 
cot a 


esc 2 a| 


cot a cot a 


= y. , cs . c . . a f v s 

cot a 


But, from the equation for the ellipse: 


y/cot a = x. 


so that 


dy _ - x 3 

da A 2 sin 2 a 


x 3 = y 3 tan 3 a. 


so that 


dy _ - y 3 tan 3 g 
da A 2 sin 2 a 


y 3 = - A 3 cot 3 a (I + A 2 cot 2 a)" 3 / 2 . 


so that 


dy _ A 3 cot 3 g (1 + A 2 cot 2 g)~ 3 / 2 tan 3 a 
da A 2 sin 2 a 


sin 2 a (1 + A 2 cot 2 a) : 


Multiplying numerator and denominator by sin a and using equation (30), 
the final expression becomes: 


N 0 A sin a (sin 2 a + A 2 cos 2 a) -3 / 2 . 


B. Parabola 

Consider next a section of a parabola with its opening, two 
units wide, along the y-axis. See Figure 5. The equation for this 
parabola is: 


y 2 = 4Px + 1, 
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FIGURE 5 


with the region of interest being for x < 0 and |y| <1. Note that 
y = ±1 when x = 0 and y = 0 when x = - ^p.. As in the case of the 
ellipse: 

y = x cot a. (37) 

Then, from equation (36): 


x 2 cot 2 a, - 4Px -1 = 0, 


so that 


4P ± (16P 2 + 4 cot 2 a ) 1 / 2 


tan 2 a : 


2P ± (4P 2 + cot 2 a ) 1 / 2 , 


but, since x < 0, tan 2 a > 0, and 2P > 0 for all a in the region of 
interest, the positive square root is discarded, and: 


t 


X 


tan 2 a 


2P - (4P 2 + cot 2 a) 1 / 2 


(38) 
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so that 


y = tan a 


2P - (4P 2 + cot 2 a) 1 / 2 


Then 


&L = 
da 


tan a 


- 4p 2 + cot 2 a) -1 / 2 (2 cot a) (- esc 2 a) 


+ sec 2 a 


2P - (4P 2 + cot 2 a) 1 / 2 


dy 1 2P - (4P 2 + cot 2 a) 1 / 2 

da sin 2 a (4P 2 + cot 2 a) 1 / 2 cos 2 a ” 


_ cos 2 a + 2P sin 2 a (4P 2 + cot 2 a) 2 - sin 2 a (4P 2 + cot 2 a) 


cos 2 a 


sin 2 a (4P 2 + cot 2 a) 1 / 2 


iz = . -2P 

da cos 2 a 


1 - 


2P 


(4P 


y + cot y 


a) T 7' ? 


jj — = dN d^ _ V 

a ~ da dy da cos 2 a 


1 - 


2P 

(4p : + cot 2 a) 1 / 2 


(39) 

(40) 


All these expressions are valid throughout the region of 
interest, except at a = 0, jt and jt/2. At a = 0, jt: cot a = “, so 

that x and y are not determined from the equations given above. 
However, these are just the points at which x = 0 and y = ±1 . At 
a = it/ 2 the situation is slightly more complicated, since Na -* 0/0 

as a jt/2. Thus to evaluate Nq, at this point, l'Hopital's Rule may 
be used as follows. Rewriting equation (39): 


dy _ 2P 
da 


(4P 2 + cot 2 a) 1 / 2 - 2P 
cos 2 a (4p 2 + cot 2 a) 1 / 2 * 
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Now let 


2P 


(4P 2 + cot 2 a)% - 2P 


= f(a) 


cos 2 a (4P 2 + cot 2 a)^ ei g(a) 


and 


(4p 2 + cot 2 a)^ = Q(a). 


Then, 


&Z_ = f(a) = 2P (Q - 2P) 
da _ g(a) q cos^ a 


Also 


Q 1 (a) = I (4p 2 + cot 2 a)~% (2 cot a) (- esc 2 a) = ~ 


and 


_ Q sin 2 q (esc 2 q) + cot Q (Q • 2 sin a cos a + Q' sin 2 a) 
Q (a) " sin 4 a 

Q 4- cos a (2Q cos a + Q 1 sin a) 

~~ sin 4 a 

so that 


Q(n/2) = 2P 
Q ' (Jt/2) = 0 

»"<*/» - If ■ i? • 
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Now 


dy _ lim f (a) _ lim f 1 (Ct) 

dal a -> jc/2 g(a) a -> n/2 g'(a) 


by 1' Hopital's Rule. 


f'(a) = 2PQ ' 


0 ! — ■» rt /2 


0 


since Q'(rt/2) = 0 and P < oo. 


g'(a) = Q' cos 2 a + 2 cos a (- sin 0!) Q — 

a 


since Q'(rt/2) = 0 and cos rt/2 = 0. Hence, 


lim f 1 (0!) _ 0 

a -> tc/ 2 g' (a) 0 


and l’Hopital's Rule may be used again. 


f"(a) = 2PQ" > — = 1 

a -4 rt/2 ** 


g"(a) = - 2Q ' cos a sin a + Q" cos 2 a - 2Q (cos 2 a - 

- 2Q ' cos a sin a 0 + 0 - 4P(0 -1) - 0 = 4P 

a -> ^2 

so that 


and 


da 



N 


a 


= % 

8P 


4P 


for a = rt/2. 


> 0 

rt /2 


sin 2 a) 


(41) 
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SECTION IV. DISTRIBUTION OF MOLECULES ACROSS THE 
Y-AXIS UPON EXIT FROM ARRAY 

A. Circular Array 

1. General Expression 

From Section II, page 4, the number of molecules per 
radian reflected from a point P^ at an angle cp with the normal is: 

N = i N cos Cp. 

a, cp 7 a y 

Now the total number of molecules per radian from P a that passes through 
a section of the opening from the origin (0, 0) to point (0, z) is 
(see Fig. 6): 



dcp 


8 



o 


\ \ cos Cp dcp = j N a 


sin S 



(o,i; 


(0,z) 


( 0 , 0 ) 


( 0 ,- 1 ) 


FIGURE 6 



and the total number that passes through this section from all points 
on the surface is: 


N 

z 



Tt 



O 


?“a 


sin 5 da. 


(42 


Thus, to calculate N 2 , the initial distribution over the surface, N a> 
and sin 5 must be expressed in terms of a. 

In the case of the circle of unit radius, the calculation of 
sin b is especially simple. From Figure 6: 


sin 5 _ sin (it - a) 
z b 


or 


sin 5 


z sin a 
b 


Also 


b = 


1 + z c 


2z cos (jt - a) 


% 


> 


and 


cos (jt - a) = - cos a 


so that 


b = (1 + z 2 + 2z cos a) 


Then 


z sin a 

( 1 + z k H- 2z cos a)% 9 


sin & 
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so that the expression for the distribution, equation (42), becomes: 


N 


Ti 

-t / 


N sin OC dOC 
CL 


(1 + z 2 + 2z cos Ct) 2 


(43) 


2. Incident Distribution Uniform over the Surface 

In this case, Nq; = ^ = constant, so that: 

at 


N 


V r sin a 

2n J (1 + z 2 + 


da 


2z cos a)' 


n z r 

o 
2n 


(1 + z 2 + 2z cos a)% 


N 

o 

2it 


(1 + z 2 - 2z)^ - (1 + z 2 + 2z) 


% 


N 

o 

2it 


“ 

± (1 - z) - ±(1 + z) 


o i 

But (1 + z + 2z cos OC ) 2 = b, which is a distance, and therefore 
positive. Since 0 < z < 1, only the positive square roots are meaning- 
ful , so that: 


N ■ 
z 


N z 
o 


Thus the exit distribution is uniform across the y-axis, for molecules 
exiting the array after the first collision. 
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3. Incident Distribution Uniform across the Y-Axis 
In this case, from equation (21), page : 

N= ^ N sin a (21) 

u 2 o 

so that, from equation (43): 


N 


z 



Nq, sin a da 

(1 + z 2 + 2z cos a)^ 



, N sm 2 a da 

1^ _o 

(1 + z 2 4- 2z cos a) 2 


N 

z 


T N z 
4 o 


it 



o 


sin 2 a da 

(1 + z 2 + 2z cos a)^ 


(44) 


This integral cannot be evaluated in closed form, but 
has been numerically integrated by a General Electric 225 Computer, 

The results of this calculation are given in Section VI of this report, 
for N c = 10,000 and several values of z. 


B. Elliptical Array 

Consider now a beam of molecules incident from the right onto 
a semi-elliptical surface as shown in Figure 7. The calculation of the 
exit distribution of those molecules which exit the array after one 
collision will depend upon the relations of the variables y, y G and z, 
as follows. 

Consider first the case shown in Figure 7, where 
y Q — 0 < 2 ^ y» I n this case the number of molecules per radian from 
Pq that passes through a section of the opening from (0, 0) to (0, z) 
is given by:- 


€ 

N z a = \ J N a cos °P d( P = \ N a ( sin € ” sin 8 ) 

& 
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where Nq; is the incident distribution function defined in Section III, 
page 43. The total number of molecules that pass through this section 
from all points on the semi-elliptical surface would then be: 


N 


z 


1 

2 


n 

/■ 


% (sin € 


sin &) da. 


(46) 


however, care must be taken to consider several different cases, and 
to define the angles 5 and e more precisely. For example, in Figure 8, 
e is less than 5 so that equation (45) would yield a negative result 
for N z the number of molecules passing through the section from 
(0, 0) to (0, z). This of course must be avoided. 

Another ambiguity is introduced when sin 6 and sin € are 
expressed in terms of a. From Figure 7: 


sin .€ _ sin 8 

z - y Q " b 



FIGURE 7 



in the case shown in Figure 9 


FIGURE 8 


FIGURE 9 


Now, a careful examination of Figures 7, 8, and 9 shows 
that these three cases are the only ones that need be considered, as 
follows. 


In the calculation of sin e and sin 8, only the triangles 
containing these angles are important,* and since the variable y is not 
a dimension of any of these triangles, its relationship to y 0 , z and 
the origin (0, 0) need not be considered. This reduces the number of 
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possible cases to six , the number of permutations of three items 
(y , z, and 0). Furthermore, half of the remaining cases are for 
negative z. Since all the arrays considered in this study are sym- 
metric with respect to the x-axis, all the results for negative z are 
identical to those for positive z, i.e. , N_ z = N z . It therefore remains 
to consider only the three cases shown in Figures 7, 8, and 9. 

Case 1. 


y Q — 0 < z (See Figure 7.) 

N z ,a = I N a (sine ’ sin&) 


where 


sin € 


(z - y n ) sin 9 
b 


and 


sin 


5= - yn sin 9 

S 


so that 


N z,a = I N a sin 9 


jz - y 0 ) , 2a 

b s 


(47) 


Case 2. 


0 < z < y c (See Figure 8.) 


N z,a 2 N a 


(sin 5 - sine) 


where now 


sin 5 


y Q sin (-9) _ - yq sin 6 
s 


s 



so that 


K z,a * 2 N a sln 0 


-ZS. . (z - Vn) 
s b 


which is the negative of the result for Case 1 above. 
Case 3. 

0 < y Q < z (See Figure 9.) 

N z,a = I N a (sin s + sin e) 


(48) 


where now 


sin 5 


y Q sin (-9) = - y Q sin 9 

s s 


and 


sin € 


(z - Yn) sin (jt - (-9)) = (y 0 - z) sin 9 
b b 


so that 


N z,a= 2 N a sin 9 


iza + (y ? .,.-,. 2 ! 

s b 


(49) 


which is the same as the result for Case 2. 

To take care of the difference in sign between Cases 1, 2, and 
3, it is necessary only to insert absolute value signs. Thus, the 
desired expression for the exit distribution is: 


59 



o 


(50) 


where the quantities Nq, y G , b, s and 0 must now be expressed in terms 
of OC (z is the independent .variable) . The incident distribution, Nq, 
for an ellipse, has been calculated in Section III-A, page 43, for an 
incident flux uniform across the y-axis. It therefore remains only to 
calculate yo, b, s, and 0: 

From Figure 7: 


tan 9 = - x/(y - y Q ) 


and from Figure 8 or 9: 

tan (-0) = - x/(y Q - y) , 


so that, in all cases: 


y D = y + x cot 0. 


(51) 


Now the equation for the ellipse is: 


x , 2 - 

j5 + y - i 


SO uiaL 


2 x dx 
A 2 


+ 2y dy = 0, 


and 


dy _ - x 
dx yA 2 ’ 

but 


- tan e, 
dx 



so that 


tan 0 = - x/yA 2 


Also from Figure 7: 


but 


tan (a - j) = y/-x, 


tan (a - -|) 


sin (a - “) 
cos (a - j) 


- cos a 
sin a 


so that 


cot a = y/x. 


From Figure 8 or 9: 


tan OL = -x/-y, 

so that in all cases: 


and 


x = y tan 
y = x cot 



Then from equation (52): 


and 


tan 0 = - tan a/ A 2 

0 = tan " 1 (- tan a/ A 2 ) 


cot a. 
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so that, from equations (53) and (51): 


y Q = y + y tan a cot 9 




y *< 1 + tan a • cot 


= y -|"l + tan a 


tan -1 (- tan a/ A 2 ) 




tan a 


y D = y (1 - A 2 ). 


The expression for y is given in Section III-A, page 45: 


y = - A cota (1 + A 2 cot 2 CC)~%. 


Now, from Figure 7: 


b = + <s/ (y - z) fc + x c 


v 2 , „2' 


and from Figure 8 or 9: 


b = + «/"( z - y) 


2 2 * 
+ X 


Also, from Figure 7, 8, or 9: 


s = + + y 


.2 . „ 2 ' 


so that, inserting equations (55), (57) and (58) into equat 


' \ I H a 


- Z - ^-- - A 2 ) ) sin 0 


7y - z) y + x 2 ' n/^FTT 51 


(55) 

(56) 

(57) 

(58) 

ion (50): 


da 


(59) 



where, from page 43: 


N a “ \ N o A sin a ( sin2 os + A 2 cos 2 a)" 3 / 2 ; 
from equation (56): 


y - - A cot a (1 + A 2 cot 2 a)"^ ; 
from equation (53): 

x = y tan a; 

and, from equation (54): 

0 = tan" 1 (- tan a/A 2 ) . 


z is of course the independent variable; A is the constant of the 
ellipse; N 0 is the number of incident molecules. 

Now, at a = 0 and a = jt, ,the expression for y (and therefore 
x) becomes indeterminate, since cot 0 = cot n = oo . However, from the 
way in which the ellipses are defined, at a = 0, y - -1 and x = 0, while 
at a = it, y = +1 and x = 0. Also, at a = rt/2, tan a = «>, but this simply 
means that 0 = ir/2 also (which can be seen from Figure 7, 8, or 9). 

Note that for A = 1, equation (59) reduces to the equation 
for the circle, as if should: 

N a = 2 sin a ( sin2 a + cos2 C0~ 3 / 2 = N q sin a 


y = - cot a (1 + cot 2 a)"^ 
>t a (i 


= - col 


cos 2 q N ^ 
-in 2 a 
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= - cot a (- in 2 OL±cosf c<Y% 

sin 2 a 


cot a (sin a) 


-X 


y = - cos a 


x = y tan a = - cos a tan a 


x = - sin a 


0 = tan -1 (- tan a) 


tan 9 = - tan a 


sin 0 - tan 9 




tan 2 9 


2 


tan a (1 + tan^ a) 


tan a 1 + 


sin 2 aV % 
cos 2 a 


- tan a (± cos a) 


sin 0 * i sin a, 


but since only the absolute value enters into the calculation: 


sin 0 = sin OC 


so that: 


it 


N = ~r 


2 N 0 sin a 


- z sin a 


/(- cos a - sin^ a’ 


da 
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-N* I 


sin 2 g da 




cos 2 a + 2z cos a ~+ z 2 + _ sin 2 "a l 


7t 


1 „ T sin 2 a 

= 4 N o Z 7 7“ 

H J (1 + z 2 + 


da 


(1 + z + 2z cos a) 


which is identical to equation (44), page 54. 

Equation (59) for the distribution of molecules exiting 
a semi-elliptical array after one collision is not integrable in closed 
form, but has been numerically integrated by the GE 225 Computer, for 
an incident flux uniform across the y-axis*. The results are given in 
Section VI of this report for N 0 = 10,000 and several values of A and z. 

C. Parabolic Array 

By referring to Figure 7 of the previous section, it can be 
seen that the expression for the exit distribution from a parabolic 
section will be very similar to that for the semi-elliptical array, 
except that the expressions for x, y, 0, y Q and Nq; will be different, 
as follows. 

The equation for the parabola is: 
y 2 = 4Px + 1, 


and from Section III-B, page 46: 


' y = x cot a 


( 60 ) 


so that 


x 2 cot 2 a - 4Px -1=0 
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and 


x 


4P ± (16P 2 + 4 cot 2 q)^ 
2 cot 2 a 


2P - (4P 2 + cot 2 q)% 
cot 2 a 


(61) 


the negative sign being chosen because the domain of interest is for 
x < 0 and P > 0. Again using the equation for the parabola: 


y 2 = 4Px + 1 

2ydy = 4Pdx 

dy 2P 
tan 0 = —• = — 
dx y 

0 = tan" 1 (2P/y). (62) 


Now, 


y Q = y + x cot 0 


(63) 


as before, and the. expression for N a (for incident beam uniform across 
the y-axis) is given by equation (40) page 48. The expression for the 
exit distribution is then given by 


N z 



N P 

o 

cos 2 OC 


2P 

(4P 2 + cot 2 a)^ 


y o - z 

n/ ( y - z) 2 + x 2 ' 



| s in 0 | da 


(64) 
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where x, y, y Q and 0 are given by equations (61), (60), (63) and (62) 
above, z, of course, is the independent variable; P is the constant 
of the parabola; N 0 is the number of incident molecules. As in the 
case of the ellipse, care must be taken at Oi ~ 0, jc/ 2, jt. At a = 0, jt 

cot a = oo, but again x « 0 and y = ±1 at these points. At CL = jt/2, it 

has already been shown in Section III-B, page that Nqj = N 0 /8P. 

Again the expression for N z is not integrable in closed form, 

but results of the numerical integration are given in Section VI of 
this report. 


SECTION V. DISTRIBUTION OF MOLECULES ALONG THE 
CENTER LINE (X-AXIS) OF AN ELLIPSE AND A PARABOLA 

A. General Expressions 

The distribution of molecules along the center line of an 
ellipse (or parabola or hyperbola) can be found by calculating the 
number of molecules which passes through a section of the x-axis 
between x x and X£ (Fig. 10). 



FIGURE 10 
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Using the same argument as in Section IV (page 51), this number would 
be given by: 



N a (sin € - sin 5) da 


(65) 


if now the angles € and 5 are defined in Figure 10. Since all the arrays 
considered are symmetric with respect to the x-axis, the integration 
extends from ft/2 to ft and is doubled to obtain the total contribution. 

T}xe problem now is to express sin € and sin 8 in terms of the 
independent variables a and x x . For purposes of numerical evaluation 
of the integral, x 2 has been taken to be x x 4- A/8 for the ellipse and 
xi + 1/48P for the parabola., The initial distribution function, Na, 
has been calculated for the ellipse and parabola in Section III above. 

Now, from Figure 10, it is clear that there will again be a 
number of cases to consider, depending upon the relative positions of 
x, Xq, x x and xg. (The position of the origin does not matter, since 
none of the pertinent triangles contain the origin.) If all permutations 
of these variables were to be considered, there would be 4J = 24 possible 
cases. Fortunately, this number can be reduced, as follows: (a) as 

mentioned above, x 2 is always greater than x x . This eliminates half of 
the 24 cases* (b) An examination of ellipses and parabolas shows that 
x < Xq in all cases. Since this condition is independent of (a), 6 of 
the remaining 12 cases are eliminated. There are, therefore, six cases 
to consider: 


1. x < Xi < x 2 < x c 

2. X < X X < Xq < x 2 

3. x < x Q < x x < x 2 

4. x x < x 2 > x < x Q 

5. Xx < X < X 2 < Xq 

6. x x < x < Xq < x 2 . 
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Also 


I 


sin 6 sin (I ~ e ) 
c r 

The desired quantity is 


(sin e - sin B) 




e. 


Expressions must now be obtained for a, c, d, r and 0 in terms of the 
independent variables a, x x and x 2 . From the definitions of a, c, d 
and r (see above), this means that x 0 , x and y (and 0) must be expressed 
in terms of a, x x and x 2 . But this has already been done for x, y and 0 
in Section IV for the ellipses and parabolas. The remaining quantity, 
x 0 , will be calculated later (see page 72). 

Case 2. x < x x < x 0 < x 2 (See Figure 12). 



FIGURE 12 
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From Figure 12: 


sin £ 
a 



or 


Also 


or 


sin e 


a cos 9 
d 


in a sln (f - e ) 


sin 8 


c cos 6 
r 


Now, from Figure 12, it is clear that the desired quantity in this 
case is not (sin e - sin 8), but is (sin e + sin 8). This may be 
obtained by using the same expression as in Case 1 above for 
(sin e - sin 8), namely: 



the negative sign being taken care of, since c is negative in this case. 
Note that all angles are considered positive, regardless of their 
direction. 

Case 3. x < xq < x x < x 2 (See Figure 13). 


in « sln (j + 9 ) 


or 


a cos 0 


sin € - 


d 


FIGURE 13 


In this case, the desired quantity is (sin & - sin e) . 
Again, if the original (Case 1) expression is used, the correct answer 
is obtained, i.e. , 




cos 0 = (sin 5 - sin e) 


Case 4. x x < x 2 < x < x Q . 


At this point it is clear that the relative position of x 
does not enter into the calculation of (sin e - sin 6). Case 4 is 
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cos 9 
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and 


s = (x 2 + y 2 ) l/s . 


so that 


= _ (x g + y 2 ) 1 / 2 sin (q + e) . 

° cos 0 > 

Equation (65), with the appropriate substitutions from 
equations ( 66 ), (67) and ( 68 ) and the expressions for Nq;, . 0 * x and y 
from Section III, yields the distribution N x along the x-axis as a 
function of xi. Calculations for specific cases follow below. 

B. Ellipse 

From page 44, equations (32), (33) and (34): 
x = - A (1 + A 2 cot 2 a )’ 1 / 2 , 


and 

y = x cot a. 

From page 60, equation (54): 

0 = tan -1 (- tan a/ A 2 ). 
From page 46, equation (35): 


% = 7 N 0 A sin a(sin 2 a + A 2 cos 2 a )" 3 / 2 
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where N Q is an arbitrary number of incident molecules and A is the 
"shape parameter" of the ellipse: 


x 

A 2 


+ T 


1 . 


Summarizing, then: (x 1} x 2 , A, a, N 0 are independent variables), 


e = tan- 1 I - •Sf§. a 


x = - A (1 + A 2 cot 2 a)” 1 / 2 


y = x cot a 

Nq; = N 0 A sin a(sin 2 a + A 2 cos 2 cc)“ 3 ^ 2 
x Q = - ( x2 + y 2 ) l/2 sin (a + 6)/cos 0 


a = x 0 - Xi 

b = ly 2 + (x Q - x) 2 ] 1 / 2 
c = x 0 - x 2 

d = [y 2 + (xx - x) 2 ] 1 / 2 
r = [y 2 + (x 2 - x) 2 ] 1 / 2 
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C. Parabola 

From page 47 , equations (38) and (37): 
x = tan 2 a [2P - (4P 2 + cot 2 a) 1 / 2 ] 
y = x cot a. 

From page 48, equation (40): 

N P 

% = ° g ■ [1 - 2P (4P 2 + cot 2 q;)” 1 ' 2 ] 

^ cos^ a • 

where P is the shape parameter of the parabola 
y 2 = 4Px + 1 . 

Summarizing: 

x = tan 2 a [2P - (4P 2 + cot 2 a) 1 / 2 ] 

y = x cot a 

8 = tan -1 (2P/y) 

N P 

R-y = -2 -b — [1 - 2P (4P 2 + cot 2 a)" 1 / 2 ]. 
a cos a 

The remaining steps are identical to those for the ellipse. 
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SECTION VI. NUMERICAL RESULTS 

A. Fraction of Molecules Making Collision with a Semi- 

Circular Surface. 

From the expressions derived for distributions over the semi- 
circular surface (page 39 ), the total number of molecules making the 
n th collision (with the surface) can be found by integrating N^ n ) from 
a = 0 to it. It is easier, however, to calculate 2 x N^ n ) at p ~ 0, which 
is the same quantity (see Figure 1, page 2). The results tollow: 

1. Incident Distribution Uniform over the Surface 


N 


2; x =2x “ (*/2) = N c 


2 x = 2 * T (f - 0 ' N ° 3 °' 3634 


N r 


2 x N^^ = 2 x -r 2 - 
P=0 4jt 


2 * + (- ? + £ 


[) - M 0 - 0.1359 M c 


n\ N 

2 x N ^ — 2 x — — — 
x w p =0 16n 


8it + ( - 38 + -y- 


- No )- °- 0505 B o 


2 * * 2 * 6^ 


N r / c-i 23 

32« + [ - 187 + ■— " -V + 24 


= 0.0188 N r 


2. Incident Distribution Uniform across the Y-Axis 
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2 x N 


(3) = 
p=o 


2 x N 


(4) = 

P=0 


x n (5 ) = 
x w p=o 






= 0.1271 N 0 




- ^ * °* 0471 N o 


2 x !° ((V _ 

48 ^ 10J5 72 16 



0.0175 N 0 


From these results, the fraction of molecules exiting the 
array after n collisions can be calculated: 


a 


n 


(# making 


collision) - (# making (n+l) th collision) 
(# making nth collision) 


The results of the trivial calculation are given in Table I. 

B. Distributions across Openings and Center -Lines of Arrays 
after One Collision, with Incident Distribution Uniform 
Across Y-Axis. 



FIGURE 14 • 
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In the analysis of the preceding sections, expressions 
were derived for the following: 


N = number of molecules that pass through a section of 
the opening between (0, 0) and (0, z), after one 
collision with the surface. 

N = number of molecules that pass through a section of the 
x-axis between (x l5 0) and (x x +*0.1, 0), after one 
collision with the surface. 


Programs were written for a General Electric 225 Computer 
to evaluate the integrals for these expressions, for several different 
geometries. The results of these calculations are given in Tables II 
through V, and are shown graphically in Figures 15 through 30. 

In Tables II and III (for exit distributions), the indepen- 
dent variable z is given in column (A), and N z is given in column (B) . 

But the quantity of greater interest is the number between z and z - 0.1. 
This number, Naz( 225), is given in column (C) . The numbers in column (D) , 
NAz(7090), are the corresponding values obtained by the Monte Carlo method, 
using the IBM 7090 Computer 1 . The difference between these two. A, is 
given in column (E) and the percentage difference in column (G) . The 
probable error , P. E., given in column (F) , is calculated as follows. 

In general, if N measurements of a quantity are taken 
and p is the probability of obtaining a given measurement, then the 
standard deviation of a measurement is given by: 

a =^N p (l - p)^ 


and the probable error , P. E., is given by P. E. = 0.6745a. The 
probability, p, is found from the 225 results (see next paragraph) and 
N = 10,000 is the number of "measurements." For example, the first 
value of N^z (for z = 0.1) in Table II is 957. This means that the 
probability, p, of a molecule falling into the section between z = 0 
and z = 0.1 is 957/10,000 = 0.0957 (since the incident number of 
molecules is 10,000). Then (1 - p) - 0.9042 and 


a =4 (10,000) (0.0957) (0.9042)' = n/ 865. 3194 1 = 29.42, 

1 This method and the results have been reported in Marshall 
Technical Paper MTP-AER0-62-53 by J. 0. Ballance, W. K. Roberts, and 
D. W. Tarbell, and were presented at the Cryogenic Engineering Conference, 
August 14-16, 1962. 
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so that 


P. E. = (0.6745) (29.42) = 19.8 = 20. 


In Tables IV and V (for center-line distributions), the 
independent variable, x, is given in column (A) and N x (225) in column (B) . 
Since N x is already for a fixed interval, no subtraction is necessary. 
N x (7090) is given for the same interval in column (C); the difference, 

A, in column (D) ; the % difference in (F) , and the probable error, P. E., 
in column (E) . 

In all cases, the GE 225 Program employed the Simpson 
Rule to evaluate the integral. Since the accuracy of these evaluations 
depends on the size of the increment, A3, used in the Simpson Rule, this 
increment was taken to be just small enough so that the GE 225 result 
contained as many significant figures as the 7090 result. The value 
used for Ax ranged from 1° down to 0.005°, depending on the particular 
case; in each case, the increment used was made smaller until two results 
were identical (to the same number of significant figures as there were 
in the 7090 result). The numbers given in column (C) for N Az (225) can 
therefore be considered exact, with the error in the Monte Carlo results 
only. 


SECTION VII, DISCUSSIONS AND CONCLUSIONS 
A. Surface Distributions 

The distribution of molecules on the inner surface of a semi- 
circle (or a long cylinder in three dimensions) has been derived as a 
function of the number of collisions made by the molecules. Initial 
distributions considered were for completely uniform incidence upon the 
surface and for a directed, uniform flow on the entrance to the surface. 
Figure 31 graphically presents the percentage of particles colliding per 
radian for a semicircle and the two initial distributions considered. 

It is seen that, for either initial distribution, after two collisions 
the distribution on the surface becomes approximately constant. Of 
course, the percentage of particles making the third, fourth, and fifth 
collisions becomes quite small (i.e., 1.75% make 5 collisions). 

Chahine [5] has derived the total number flux incident at any point on 
a semicircular surface per unit time and unit area with a directed uni- 
form flow. By summing the distributions for the five collisions this 
same flux may be found. This summation is within 2% of the Chahine 
result. Since approximately 99.4% of the incident particles have exited 
after 5 collisions, this difference is not significant and shows that, 
for this case, less than 2% of the flux results from particles making 
more than 5 collisions. 
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B. Exit and Center Line Distributions 


The distributions of particles exiting the surface- array and 
along the center line of symmetry of the array have been derived* Due 
to the complexity of these derivations, analytical solutions were 
obtained only for single collisions with the surface; however, computer 
programs can consider any number of collisions. Early in the study it 
was thought that . concentration of particles (focusing) would occur at 
the exit of the array. No such focusing was observed, but the distribu- 
tions at the exit and along the center line does yield very interesting 
results. For example, for an elliptical surface with the parameter A 
equal to two, it is seen from the Table II and Figure 18 that 50 % of the 
incident particles cross the inner 50% of the center line after the 
first collision. A practical application of the result is shown in 
Figure 22. Here is a cryogenic pumping array which could be used for 
highly directed flow fields such as in a low density wind tunnel, a 
rocket exhaust test chamber, or a rocket sounding probe for atmospheric 
sampling. Radiation shields reduce the radiant heat losses of the 20°K 
surface and also provides some impedance to molecules from leaving the 
array* Collection efficiency (that is, the percentage of particles 
condensed by the 20°K surface compared to the total number of particles 
entering the array) should be very high (hopefully, 80-90 percent). 

C. Monte Carlo Computer Technique g 

The applicability and accuracy of Monte Carlo computer techni- 
ques for the study of free molecular flow has been shown. These techni- 
ques are necessary for this type of study since the analytical solutions 
become quite complex. As future plans in the study are for extension of 
the program into transition flow where intermolecular collisions will be 
considered, it is felt that only through the use of Monte Carlo methods 
can meaningful studies be made. 
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APPENDIX A 
Useful Integrals 


r a u , . 2 ^ 

/ sm — cos — dCK = sin ~ 


(1) 


J 


r .,.2 a . 1 


sin 2 da = “ sin a) 


( 2 ) 


/ 


Qf 1 

cos 2 ~2 doc - ~2 + sin 


( 3 ) 


I 


a sin j cos y da = (sin a - a cos a) 


( 4 ) 


/ 


a 2 sin ^ cos ^ da = a sin a + cos a - — a 2 cos a 


( 5 ) 


sin a sin ^ cos ^ da = (a - sin 2a) 


( 6 ) 


I 


cos a sin ^ cos ^ da = sin 2 a 


( 7 ) 



_1 

2 


(i a 2 +' a sin a + cos a) 


( 8 ) 



APPENDIX A (Cont'd) 


a 2 


COS 


p GC - 

tt da = 




sin Ot + OL cos a 


(9) 


sin a cos 2 



4 a 
cos* — 


( 10 ) 


cos a cos 2 J da = ^ (^a + sin a + sin a cos 2 ^ + X sin 2a 


2 4 


(ID 


Q 1 / | 

a sin 2 — da = j ( — — a sin a - cos a 


( 12 ) 


a 2 sin 2 J da = - (-Y - 1 ) Sin a - a cos a 


(13) 


. 2 a J/v . 4 a 
sina sin j da = sm* 


(14) 


cos a sin 2 ^ da = j a + -| sin a + sin a cos 2 ~ sin 2a 


(15) 


a , „ . a 

cos — da = 2 sm — 


(16) 


sm — da = 


o a 

- 2 cos — 


(17) 
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APPENDIX A (Cont'd) 


QL CC CC 

CC sin — da = 4 sin ^ - 2a cos — 


(18) 


I 


a cos ^ da = 4 cos ^ + 2a sin ^ 


(19) 


I 


cos ^ da = 4a cos ^ + fa 2 - 8 ^ sin ^ 


( 20 ) 


/ 


sin ^ da = 4a sin ^ - fa 2 - 8^ cos ~ 


( 21 ) 


^ a 3 cos ^ da = (12a 2 - 96) cos + (2a 3 - 48a) sin — | 


( 22 ) 


/ 


a 3 sin ^ da = (12a 2 


96) sin | - (2a 3 - 48a) cos | 


(23) 



TABLE I 

FRACTION OF MOLECULES THAT EXIT 
SEMICIRCULAR ARRAY AFTER THE n th COLLISION 




A 

B 


n 


number making 
n 11 * 1 collision 
4- N 0 

number exiting 
array after . 
n fc h collision 
4-N 0 

ii 

>|W 

i 


1.0000 

0.6366 

0.6366 

2 


0.3634 

0.2275 

0.6260 

3 


0.1359 

0.0854 

0.6284 

4 


0.0505 

0.0317 

0.6277 

5 


0.0188 

- 

- 

1 


1.0000 

0.6667 

0.6667 

2 


0.3333 

0.2062 

0.6187 

3 


0.1271 

0.0800 

0.6294 

4 


0.0471 

0.029B 

0.6285 

5 


0.0175 

- 

- 
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TABLE II 

ELLIPSE EXIT DISTRIBUTION 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) 

(G) 

Z 

N z 

N Az (225) 

N^(7090) 

A 

P.E. 

A(%) 

A = 0.50 







0.1 

956.5 

957 

936 

21 

20 

2.2 

0.2 

1909.7 

953 

924 

29 

20 

3.1 

0.3 

2856.1 

946 

961 

15 

20 

1.6 

0.4 

3791.3 

935 

934 

1 

20 

0.1 

0.5 

4710.1 

919 

924 

5 

19 

0.5 

0.6 

5604.9 

895 

893 

2 

19 

0.2 

0.7 

6464.7 

860 

935 

25 

19 

2.7 

0.8 

7271.1 

806 

827 

21 

18 

2.5 

0.9 

7989.5 

718 

691 

27 

17 

3.9 

1.0 

8527.5 

538 

529 

9 

15 

1.7 

A = 1.00 







0.1 

784.4 

784 

805 

21 

18 

2.7 

0.2 

1562.9 

779 

729 

50 

18 

6.4 

0.3 

2329.4 

767 

7.91 

24 

18 

3.1 

0.4 

3077.4 

748 

797 

49 

18 

6.6 

0.5 

3800.1 

723 

744 

21 

17 

2.9 

0.6 

4489.5 

689 

640 

49 

17 

7.1 

0.7 

5136.3 

647 

624 

23 

17 

3.6 

0.8 

5728.8 

593 

568 

25 

16 

4.2 

0.9 

6249.9 

521 

477 

44 

15 

8.4 

1.0 

6666.7 

417 

433 

16 

13 

3.8 
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TABLE II (Cont'd) 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) 

(G) 

Z 

N z 

N Az< 225) 

O 

o\ 

O 

r^. 

'w' 

S3 

A 

P.E. 

A(%) 

A = 1.50 







0.1 

600.7 

601 

597 

4 

16 

0.7 

0.2 

1196.7 

596 

583 

13 

16 

2.2 

0.3 

1783.1 

586 

583 

3 

16 

0.5 

0.4 

2355.2 

572 

676 

104 

16 

15.4 

0.5 

2907.8 

553 

554 

1 

15 

0.2 

0.6 

3435.8 

528 

557 

29 

15 

5.2 

0.7 

3933.2 

497 

500 

3 

15 

0.6 

0.8 

4393.2 

460 

483 

23 

14 

4.8 

0.9 

4806.7 

414 

421 

7 

13 

1.7 

1.0 

5157.6 

351 

348 

3 

12 

0.9 

A = 2.00 







0.1 

460.5 

461 

486 

25 

14 

5.4 

0.2 

917.8 

457 

468 

11 

14 

2.4 

0.3 

1368.6 

451 

452 

1 

14 

0.2 

0.4 

1809.7 

441 

418 

23 

14 

5.2 

0.5 

2237.7 

428 

415 

13 

14 

3.0 

0.6 

2649.1 

411 

385 

26 

13 

6.3 

0.7 

3040.1 

391 

375 

16 

13 

4.1 

0.8 

3406.3 

366 

361 

5 

12 

1.4 

0.9 

3742.1 

336 

334 

2 

12 

0.6 

1.0 

4037.7 

296 

288 

8 

11 

2.7 
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TABLE III 

PARABOLA EXIT DISTRIBUTION 


(A) 

(B) 

(c) 

(D) 

(E) 

(F) (G) 

z 

N 

z 

1^(225) 

N^(7090) 

A 

P.E. A(7.) 


P = 0.5000 


0.1 

942.1 

942 

968 

26 

20 

2.8 

0.2 

1879.6 

938 

930 

8 

20 

0.9 

0.3 

2808.1 

929 

941 

12 

20 

1.3 

0.4 

3723.5 

915 

889 

26 

19 

2.8 

0.5 

4621.5 

898 

927 

29 

19 

3.2 

0.6 

5498.7 

877 

866 

11 

19 

1.3 

0.7 

6351.4 

853 

816 

37 

19 

4.3 

0.8 

7176.4 

825 

848. 

23 

19 

2.8 

0.9 

7970.2 

794 

766 

28 

18 

3.5 

1.0 

8727.6 

757 

771 

14 

18 

1.8 


P = 0.2500 


0.1 

748.5 

749 

732 

17 

18 

2.3 

0.2 

1492.1 

744 

721 

23 

18 

3.1 

0.3 

2225.9 

734 

734 

0 

18 

0 

0.4 

2945.4 

720 

715 

5 

17 

0.7 

0.5 

3646 . 6 

701 

701 

0 

17 

0 

0.6 

4325.7 

679 

678 

1 

17 

0.1 

0.7 

4979.4 

654 

645 

9 

17 

1.4 

0.8 

5604.7 

625 

608 

17 

16 

2.7 

0.9 

6198.7 

594 

611 

17 

16 

2.9 

1.0 

6756.8 

558 

547 

11 

15 

2.0 
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TABLE III (Cont'd) 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) (G) 

Z 

N 

z 


/-N 

o 

VO 

o 

A 

P.E. A(%) 


P = 0.1250 


0.1 

439.4 

439 

453 

14 

14 

3.2 

0.2 

877.0 

438 

422 

16 

14 

3.7 

0.3 

1310.9 

434 

403 

31 

14 

7.1 

6.4 

1739.5 

429 

428 

1 

14 

0.2 

0.5 

2161.1 

422 

394 

28 

14 

6.6 

0.6 

2574.0 

413 

392 

21 

13 

5.1 

0.7 

2976.7 

403 

415 

12 

13 

3.0 

0.8 

3367.6 

391 

393 

2 

13 

0.5 

0.9 

3745.0 

377 

420 

43 

13 

11.4 

1.0 

4106.6 

362 

400 

38 

12 

10.5 


P = 0.0833 


0.1 

288.3 

288 

283 

5 

11 

1.7 

0.2 

575.9 

288 

279 

9 

11 

3.1 

0.3 

862.0 

286 

275 

11 

11 

3.8 

0.4 

1146.0 

284 

264 

20 

11 

7.0 

0.5 

1427.0 

281 

265 

16 

11 

5.7 

0.6 

1704.5 

278 

258 

20 

11 

7.2 

0.7 

1977.6 

273 

260 

13 

11 

4.8 

0.8 

2245.6 

268 

291 

23 

11 

8.6 

0.9 

2507.6 

262 

256 

6 

11 

2.3 

1.0 

2762.3 

255 

272 

17 

11 

6.7 
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TABLE IV 

ELLIPSE CENTER LINE DISTRIBUTION 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) 

X 

N x (225) 

N x (7090) 

A 

P.E. 

A(%) 

A = 0.5000 

-0.5000 

405 

403 

2 

13 

0.5 

-0.4375 

385 

373 

12 

13 

3.1 

-0.3750 

364 

393 

29 

13 

8.0 

-0.3125 

343 

358 

15 

12 

4.4 

-0.2500 

322 

301 

21 

12 

6.5 

-0.1875 

301 

266 

35 

12 

11.6 

-0.1250 

280 

265 

15 

11 

5. '4 

-0.0625 

260 

278 

18 

11 

6.9 

A = 1.0000 
- 1.000 

872 

866 

21 

19 

2.5 

-0.875 

785 

793 

8 

18 

1.0 

-0.750 

698 

677 

21 

17 

2.4 

-0.625 

616 

640 

24 

16 

2.6 

-0.500 

538 

566 

28 

15 

2.8 

-0.375 

465 

438 

27 

14 

3.0 

-0.250 

399 

389 

10 

13 

3.3 

-0.125 

340 

336 

4 

12 

3.5 
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TABLE IV (Cont'd) 


(A) 

• (B) 

(C) 

(D) 

(E) 

(F) 

X 

N x (225) 

N (7090) 

X 

A 

P.E. 

A(%) 

A = 1.5000 

-1.5000 

1318 

1299 

65 

22 

1.8 

-1.3125 

1109 

1082 

27 

21 

1.9 

-1.1250 

928 

959 

31 

20 

2.2 

-0.9375 

771 

738 

33 

18 

2.3 

-0.7500 

634 

618 

16 

16 

2.5 

-0.5625 

515 

541 

26 

15 

2.9 

-0.3750 

412 

394 

18 

13 

3.2 

-0.1875 

326 

333 

7 

12 

3.7 

A = 2.0000 

-2.00 

1719 

1695 

17 

25 

1.5 

1.75 

1345 

1370 

25 

23 

1.7 

1.50 

1067 

1002 

65 

21 

2.0 

1.25 

847 

866 

19 

19 

2.2 

1.00 

668 

676 

8 

17 

2.5 

0.75 

518 

525 

7 

15 

2.9 

0.50 

393 

372 

21 

13 

3.3 

0.25 

293 

314 

21 

11 

3.8 



TABLE V 


PARABOLA CENTER LINE DISTRIBUTION 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) 

X 

N x (225) 

N x (7090) 

A 

P.E. 

A(%) 


P m 0,5000 


-0.50000 

284 

260 

24 

11 

8.5 

-0.45833 

276 

243 

33 

11 

12.0 

-0.41666 

267 

298 

31 

11 

11.6 

-0.37500 

257 

250 

7 

11 

2.7 

-0.33333 

248 

261 

13 

10 

5.2 

-0.29166 

238 

219 

19 

10 

8.0 

-0.25000 

229 

268 

39 

10 

7.0 

-0.20833 

220 

• 214 

6 

10 

2.7 

-0.16666 

210 

212 

2 

10 

1.0 

-0.12500 

201 

198 

3 

9 

1.5 

-0.08333 

193 

179 

14 

9 

7.3 

-0.04166 

184 

188 

4 

9 

2.2 


P = 0.250 


- 1.0000 

615 

589 

26 

16 

4.2 

-0.9166 

568 

571 

3 

16 

0.5 

-0.8333 

525 

549 

24 

15 

4.6 

-0.7500 

487 

489 

2 

15 

0.4 

-0.6666 

453 

477 

24 

14 

5.3 

-0.5833 

421 

408 

13 

14 

3.1 

-0.5000 

391 

409 

18 

13 

4.6 

-0.4166 

362 

334 

28 

13 

7.7 

-0.3333 

335 

346 

11 

12 

3.3 

-0.2500 

308 

303 

5 

12 

1.6 

-0.1666 

282 

261 

21 

11 

7.4 

-0.0833 

258 

238 

20 

11 

7.8 
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TABLE V (Cont'd) 


(A) 

(B) 

(C) 

(D) 

(E) 

(F) 

X 

N x (225) 

N x (7090) 

A 

P.E. 

A(%) 

P = 0.1250 

-2.000 

1169 

1116 

53 

22 

4.5 

-1.833 

941 

924 

17 

20 

1.8 

-1.666 

798 

812 

14 

18 

1.8 

-1.500 

699 

724 

25 

17 

3.6 

-1.333 

623 

642 

19 

16 

3.0 

-1.166 

562 

615 

53 

16 

9.4 

- 1.000 

509 

475 

34 

15 

6.7 

-0.833 

461 

458 

3 

14 

0.7 

-0.666 

414 

445 

31 

13 

7.5 

-0.500 

368 

371 

3 

13 

0.8 

-0.333 

322 

326 

4 

12 

1.2 

-0.166 

276 

250 

26 

11 

9.4 


P = 0.0833 


-3.0000 

1564 

1553 

11 

25 

2.0 

-2.75 

1101 

1131 

30 

21 

2.7 

-2.50 

886 

880 

6 

19 

0.7 

-2.25 

756 

780 

24 

18 

3.2 

-2.00 

666 

654 

12 

17 

1.8 

-1.75 

598 

614 

16 

16 

2.7 

-1.50 

541 

520 

21 

15 

3.9 

-1.25 

491 

463 

28 

15 

5.7 

-1.00 

443 

454 

11 

14 

2.5 

-0.75 

392 

423 

31 

13 

7.9 

-0.50 

336 

324 

12 

12 

3.6 

-0.25 

274 

285 

11 

11 

4.0 
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